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Scattering  of  a  Plane  Electromagnetic  Wave 
From  a  Semicircular  Crack  in  a  Perfectly 
Conducting  Ground  Plane 


1.  INTRODUCTION 


The  problem  of  scattering  from,  and  coupling  through,  cracks  and  gaps  arises  when 
considering  antennas  and  scatterers  assembled  from  many  subsections  that  meet  to  form 
cracks  and  gaps. 

The  crack  problem  has  been  studied  by  many  people  including  Lord  Rayleigh,  who 
considered  the  problem  of  plane  wave  scattering  from  a  half-cylindrical  indentation  in  a 
ground  plane. 1  In  this  paper  he  actually  solved  the  half-cylindrical  excrescence  (protuberance) 
problem,  with  some  discussion  of  experimental  techniques  related  to  the  Indentation,  or 
crack,  problem.  The  related  problem  of  scattering  from  a  silt  in  an  infinitely  thin  screen  has 


Received  for  publication  31  March  1989 

1  Lord  Rayleigh  (1907)  On  the  light  dispersed  from  fine  lines  ruled  upon  reflecting  surfaces  or 
transmitted  by  very  narrow  slits.  Phil  Mag,  14:350. 
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b«en  extensively  studied,2,  3i  4  as  has  the  scattering  from  a  slit  in  a  thick  screen.5,  6>  7-  8 
Integral  equations,  experimental  methods,  and  pure  numerical  techniques,  such  as  finite- 
difference  time  domain,  have  been  used  to  study  the  scattering  from  Indentations  in  a  ground 
plane.9,  10,  11  The  related  problem  of  scattering  from  a  cavity-backed  slit  has  been  considered 
recently  using  dual-series  techniques.12 

This  report  presents  a  dual-series  eigenfunction  approach  to  scattering  of  electromagnetic 
plane  waves  from  a  half-cylindrical  (circular)  Indentation  in  a  ground  plane.  We  use  to  our 
advantage  the  coincidence  of  the  problem  geometry  and  a  constant  coordinate  surface  In 
circular  cylindrical  coordinates.  The  main  difference  between  our  analysis  and  typical 
eigenfunction  problems  lies  in  the  incomplete  orthogonality  of  the  sinusoids  over  a  half  space, 
and  the  use  of  two  separate  regions  of  the  problem. 

Usefulness  of  the  solution  is  primarily  in  the  small-crack  results,  which  are  valid  in  the 
case  where  the  radius  of  the  crack  is  much  smaller  than  the  wavelength  of  the  Incident  plane 
wave.  For  both  the  transverse  magnetic  and  transverse  electric  polarizations,  the  ka- 
dependence  can  be  removed  from  the  necessary  numerical  matrix  inversion  when  the  crack  is 
small.  The  matrix  inversion  need  only  be  performed  once  for  each  polarization  and  it  results 
in  a  constant  coefficient  for  each  case.  Since  the  small-crack  results  are  then  functions  of  ka, 
the  two  results  (TM  and  TE)  can  be  used  to  consider  cases  where  the  incident  field  is  not 
normal  to  the  axis  of  the  small  crack. 


2  Bouwkamp,  C.J.  (1954)  Diffraction  theory.  Rep  Prog  Phys,  17:35. 

3  Chow,  T.Y.,  and  Adams,  A.T.  (1977)  The  coupling  of  electromagnetic  waves  through  long 
slots.  IEEE  T-EMC.  19:65. 

4  Butler,  C.M..  Rahmat-Samil,  Y.,  and  Mlttra,  R.  (1978)  Electromagnetic  penetrations  through 
apertures  in  conducting  surfaces,  IEEE  T-EMC,  20:82. 

5  Kashyap,  S.C.,  and  Hamid.  M.A.K.  (1971)  Diffraction  characteristics  of  a  slit  in  a  thick 
conducting  screen,  IEEE  T  AP,  19(No.  4):499. 

6  Hongo,  K.,  and  Ishil,  G.  (1978)  Diffraction  of  an  electromagnetic  plane  wave  by  a  thick  slit, 
IEEE  T  AP.  26(No.  3):494. 

7  Harrington,  R.F.,  and  Auckland,  D.T.  (1980)  Electromagnetic  transmission  through  narrow 
slots  In  thick  conducting  screens,  IEEE  T  AP,  28:120. 

®Taflove,  A.,  Umashankar,  K.R..  Beker,  B..  Harfoush,  F.,  and  Yee,  K.S.  (1988)  Detailed  FD-TD 
analysis  of  electromagnetic  fields  penetrating  narrow  slots  and  lapped  Joints  in  thick 
conducting  screens,  IEEE  T  AP,  36(No.  2):247. 

9  White,  G.S.,  and  Feldman,  A.  (1981)  Diffraction  from  a  shallow  rectangular  groove,  Applied 
Optics,  20(No.  14):2585. 

10  White,  G.S.,  and  Marchiando,  J.F.  (1983)  Scattering  from  a  V-shaped  groove  in  the 
resonance  domain,  Applied  Optics,  22(No.  15):2303. 

11  Dominek,  A.K.,  Shamanski,  H.T.,  and  Wang,  N.  (1988)  Scattering  from  Three-Dimensional 
Cracks.  Ohio  State  University  (preprint). 

12  Ziolkowski,  R.W.,  and  Grant,  J.B.  (1987)  Scattering  from  cavity-backed  aperture,  IEEE  T  AP, 
35(No.  5):504. 


2 


2.  TRANSVERSE  MAGNETIC  POLARIZATION 


2.1  Problem  Geometry 

We  consider  a  Transverse  Magnetic  (TM)  plane  wave  that  is  polarized  in  the  z-dlrection  and 
makes  an  angle  <pinc  with  the  positive  x-semiaxis.  It  is  incident  upon  and  scatters  from  a  half- 
cylindrical  Indentation  or  crack  in  the  ground  plane,  as  shown  in  Figure  1 . 


Figure  1 .  Problem  Geometry  for  Incident  TM  Waves 

The  crack  is  of  radius  r  =  a  and  of  infinite  extent  in  the  z-direction.  All  angles  are 
measured  positive  in  the  counterclockwise  direction  starting  with  <p  =  0  along  the  positive  x- 
semiaxis.  Angles  are  in  the  range  0  <  $  <  2  n .  The  crack  is  thus  described  by  n  <  <l>  <  2n  and 
r  =  a.  The  surface  0  <  ^  <  n  and  r  =  a  is  referred  to  as  the  aperture.  It  is  the  region 
complementary  to  the  crack.  The  problem  will  be  formulated  so  that  boundary  conditions 
along  the  ground  plane  will  be  automatically  satisfied,  and  then  three  boundary  conditions 
will  be  applied  at  the  surface  r  =  cl  one  on  the  crack  and  two  across  the  aperture.  We  define 
the  r  >  a  to  be  the  exterior  region  and  r  <  a  to  be  the  interior  region,  and  the  fields  in  these 
regions  are  the  exterior  and  interior  fields  respectively. 


2.2  Exterior  Field 

The  incident  field  is  a  plane  wave  as  discussed  above.  We  write 


E<nc  =  elkr  co*^  *"*) 
z 


(1) 
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where  unit  amplitude  has  been  assumed  and  etof  time  dependence  suppressed.  This  plane  wave 
has  the  following  well-known  expansion  In  terms  of  cylindrical  Bessel  functions  13 


Elnc 

z 


=  X  tnJn(kr)e^~*  1 


(2) 


In  the  exterior  region  the  scattered  field  will  be  considered  to  be  made  up  of  two  parts.  The 
first  Is  the  reflected  wave  that  would  be  present  If  there  were  no  crack,  and  the  second  is  the 
deviation  from  this  caused  by  diffraction  at  the  crack.  Thus,  in  the  exterior  region  we  refer  to 
the  Incident,  reflected,  and  diffracted  waves.  We  are,  of  course,  most  Interested  In  the 
diffracted  field. 

The  reflected  wave  Is  also  a  plane  wave  and  is  well-known.  We  have  <pref  -2k-  <pinc  so  that 
the  reflected  plane  wave  can  be  written  as 


Z 


(3) 


In  terms  of  cylindrical  Bessel  functions  we  have 


i"j 


(4) 


For  the  diffracted  wave  In  the  exterior  region  (r  >  a,  0  <  <p  <  n)  we  write 


e'^  -  X  AJhl^(kr)  sin  nip  , 

n=0 


(5) 


where  An  is  an  unknown  modal  coefficient  and  we  have  expanded  this  outgoing  wave  in  terms 
of  Hankel  functions  (of  the  second  kind)  and  also  sin  nip.  Noting  that  because  sin  0  = 
sin  n/r  =  0,  the  diffracted  wave  vanishes  on  the  ground  plane  and  so  satisfies  the  boundary 
condition  there.  Calling  the  interior  Region  1  and  the  exterior  Region  2,  the  electric  fields 
become  respectively,  E1  and  E2.  If  we  write  the  total  field  In  the  exterior  region  as 


13  Morse,  P.M.,  and  Feshbach,  H.  (1953)  Methods  of  Theoretical  Physics,  McGraw-Hill  Book 
Company,  New  York. 
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(6) 


Z 


Edf, 


we  see  that  =  0  for  <p  =  0,  n  that  is,  on  the  ground  plane. 


2.3  Interior  Field 

In  the  interior  region  the  total  Ez-field  can  be  written 


D  J  (kr)ein*. 

n  nv  7 


(7) 


where  Dn  is  an  unknown  modal  coefficient.  However,  it  will  be  more  convenient  for  our 
analysis  if  we  write  this  as 


E*  =  Jn(kr)(Bn  cos  rup  +  Cn  sin  nip), 

n»0 


(8) 


where  Bn  and  Cn  are  the  modal  coefficients  that  will  be  determined  from  the  boundary 
conditions. 


2.4  Total  Exterior  Field 

We  now  wish  to  write  the  incident  and  reflected  Ez- fields  as  series  over  n  =  0  to  ~  rather 
than  over  all  n.  We  note  that  J.n(kr)  =  (-1  )nJn(kr), 


E^  =  -J0(kr)+  £  <"j(k r)e+<n^-^l 

ncO 


nsO 


(9) 


or 


E1™  =  -  JQ(kr)  +  2  X  lnJ„(.kr) cos  n(d>  -  0  /nc). 

nmO 


(10) 
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Similarly  for  the  reflected  wave  we  write 


=  JQ(kr)  -  2  £  inJn(fcr)cos  n($  +<p  <nc). 


rc=0 


(in 


Hence, 


re^=  2  lnJ  (kr)Lcos  n(<p -<p  fnc)  -  cos  n(0+0lnc)]. 


_(nc  _ 

E  +  E  _ 

z  z  n 

rc=0 


(12) 


If  we  use  the  trigonometric  Identity  cos(a  +  b)  =  cos  a  cos  b  +  sin  a  sin  b  we  can  write 


E^  +  =  4  £  inJ^(kr) sin  n0mcsin  nip. 

n=  0 


(13) 


Recalling  the  diffracted  field,  we  have 


E^  =  £  |  4<nJ](kr)sin  rnp^  +  AJi^\k r) sin  nip, 

rv sO 


(14) 


for  the  exterior  field.  In  this  expression,  An  is  the  modal  coefficient  that  we  desire. 


2.5  Magnetic  Field  Component 

The  ^-component  of  the  magnetic  field  is  also  of  interest  in  this  problem.  It  is  related  to 
&z  by 


H  =  — -  d  E  , 
*  cup  r  z 


(15) 


so  that  the  exterior  and  interior  H$  -fields  are 


6 


-ik 

<m 


£  [4<n  J'n(fcr)  sin  n<)><nc 

n*0 


+  A  H(2Y{kr) 

n  n  v  ' 


sin  n<p. 


(16) 


and 


Hl  =  777  ^  Jn(kr)  C°S  +  Cn  ^  "0) 

y  °v  n=0 


(17) 


where  J^(z)  =  dzJn(z). 


2.6  Boundary  Conditions 

The  required  boundary  conditions  for  solving  this  TM  problem  are 


E 


2 


z 


=  0  for  r  =  a  and  n<  <f><  2n, 


(18) 


z 


for  r  =  a  and  0  <  <p  <  n. 


(19) 


H.2  =  H*  for  r  =  a  and  0  <  tp  <  n. 

t>  0  (20) 


The  first  condition  enforces  zero  tangential  electric  field  on  the  surface  of  the  perfectly 
reflecting  crack  and  the  last  two  ensure  continuous  tangential  electric  and  magnetic  fields 
across  the  aperture,  or  (imaginary)  surface  that  is  complementary  to  the  crack. 

The  boundary  conditions  give  the  following  three  equations 


J  Jn(ka)Bn  cos  n01  +  £  JJka)Cn  sin  n<p[  =0, 

n*0  rusO 


(21) 
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(22) 


X  Jn{ka)Bn  cos  n,?2  +£  Jn(ka)Cn  sin  rvp2 

n=0  n=:0 


=  X  |4inJn(ka)  sin  ^2"%  A^^ka)]  sin  n<p2. 

n=  0 


X  J'rlka)Bn  COS  n02  +  X  Jn(ka)Cn  Sln  n<&2  (23) 

n=0  n=0 


=  £  ;4fn  J'^(ka)  sin  rup, 

n»0 


<nc  (2)'/,  v 

,  +  A  fi  (ka) 

2  n  n 


sin  n02, 


where  we  must  keep  in  mind  that  Eq.  (21)  is  valid  for  n<  <pi  <  2n  and  Eqs.  (22)  and  (23)  are 
valid  for  0  <  fa  <  n  . 

Now  define  the  following  for  convenience 


F  =  4i"j  (ka)  sin  n0,nr  +  A  Hi2\ka)  -  C  J  (ka). 

n  n'  '  nnv  /  nnv/ 


(24) 


G  =  B  J  (ka). 

n  n  n'  ’ 


(25) 


R  =  -  C  J  (ka), 

n  n  nv  ' 


(26) 


as  well  as  F ' 

n 


dF 

- ^  and  so  on. 

d(ka) 


The  Boundary  Condition  equations  are  then 


XG„cos  n0,  =  £  F,  sln  n*r  ,27, 

n*0  n=0 


for  n  <  <pi  <  2k  and 


X  Gn  cos  n<P2  =  X  Fn  sln  (28) 

n*0  o*0 


8 


(29) 


X  G'n  cos  n<t>2  =  £  F'n  sin  n02. 

n»0  n*0 


for  0  <  <p  2  <  n  . 

In  Eq.  (27)  we  make  the  change  of  variables  $>  =  <p  -  rt  so  that  we  have 


£  G  cos  n(<D  +  tt)  =  £  R  sin  n(4>  +  n), 

n=0  ri=0 


(30) 


for  0  <  O  <  n  .  Noting  that  cos  n(6  +  x)  =  (-  l)n  cos  nd  and  sin  n(0  +  x)  =  (-  l)n  sin  nd  we  have 


]T  (-  l)nGn  cos  n®  =  ]£  (-  l)nf?n  sin  n® , 

n=0  ru=0 


for  0  <  4>  <  n  .  Thus,  we  see  that  all  three  Boundary  Condition  equations  are  of  the  form 


£  cos  np  =  ]jT  Fn  sin  n0, 

n*0  rt=0 


where  0  <  <f>  <  x  . 


(32) 


2.7  Orthogonality 

We  now  make  use  of  the  following  orthogonality  relations  among  the  sinusoids 


£ 


sin  nA  sin  mipdip  =  5  m,  n  *  0, 

2  mn 


f 


sin  nip  cos  rmpd<p  = 


if  n-m  is  even 


-2n —  if  n-m  is  odd 


2  2 
n-m 


(33) 


(34) 
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We  multiply  Eq.  (32)  by  sin  m<p  and  Integrate  from  0  to  n. 


cos  n<f>  sin  m<p  d<t>  = 


sin  n0  sin  d<> . 


(35) 


Exchanging  the  sums  and  integrals  and  applying  the  orthogonality  relations  gives 


2  i 


f  mG"  1 


(36) 


where  Y  *  is  used  to  indicate  a  sum  over  n  from  zero  to  infinity  where  (n  -  m)  is  odd.  Writing 

4m*  ri 

the  even-odd  behavior  explicitly. 


F  =4 
m,  n 


oe 

I 

rm  1,3.5 


mG 


2  2 
m  -  n 


(37) 


oo 

I 

re.0,2,4 


mG 

n 


(38) 


Referring  to  our  previous  definitions  for  Fn .  Gn  and  Rn  we  write 


These  three  equations  can  now  be  used  to  solve  for  the  three  modal  coefficients  An  ,  Bn  and  C„  . 


2.8  Solution  for  Modal  Coefficients 

Although  the  modal  coefficients  of  Interest  In  this  problem  are  An ,  It  is  not  possible  to 
solve  for  them  directly.  Rather,  we  use  Eqs.  (39)  and  (40)  to  eliminate  An ,  and  after  some 
algebra  use  Eq.  (41)  to  eliminate  Cn  .  The  resulting  expression  may  then  be  solved  for  Bn  ,  and 
given  that,  we  merely  backsubstitute  to  get  An  . 

Solving  Eqs.  (39)  and  (40)  for  An  and  equating  them  gives  the  relation 


-4fmJ  (ka)Hl2)  (ka)  sin  m<|>,nc  +  C  J  (ka)H®  (ka) 

mv/mv/  mm'  7  m  ' 


+  (ka)  4m  £ 


Jt  ~  „2  2 

'  n  m  -  n 


(42) 


=  -4im  J'  (fca)//2,(ka)  sin  m/nc  +  C  J'  (ka)//t2,(ka) 


(2)  .  ~  B  J'(ka) 

+  H  2  (ka)  4m  V  a  n_  ' 
m  _  „2  _2 


n  m  -  n 


If  we  rearrange  and  use  the  Wronskian  relationship,  W 


-2( 

tfka 


we  can  write 


(4im  sin  m<plnc  -  C  )-2L 
'  m*nka 


(43) 


=  4m  y 
re  " 


2  2 
n  m  -  n 


J'(ka)H[2\ka)  -  J  (ka)Hl2] 

tv  m  '  7  n'  7  m 


(ka)  . 


Substituting  for  Cn  we  find 


(44) 


4lm  sin  m»“  =  4m  X  -A- 

*  n  m2-  n2  Jm(ka) 


*  C1^)4?  i -  j im»>] 

\  Zl  /  n  nm-n 


Finally,  rearranging  this  gives 


i f -v’-al {j.(fca) +  (¥)[jn(ka)Hi21  (ka)jm(fca) -  j;(ka)Hi21(ka)jm(ka)]}  (45) 

n  -  n  y  \  ^  / 


=  -2-  imJ  (ka)  sin  mtpinc 
m  m 


which  we  can  solve  for  Bn  .  Once  the  Bn  are  found,  the  Ai  can  be  found  from 


■■■{??(¥? 


-  4 ilJ(ka)  sin  £0<nc 


}  nj2\ka) 


(46) 


which  determines  the  diffracted  field  in  the  exterior  region. 


2.9  Small  Crack  Limit 


When  ka  -»  0  the  cylinder  functions  have  the  following  simple  forms 


J,(ka) 


(ka/2) 


(47) 


and 


< W)  -L  _1 1— 

l*  (ka/2) 


(48) 
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The  equations  that  we  must  solve  also  become  much  simpler,  and  we  can  separate  the  ka- 

dependence  from  the  necessary  numerical  matrix  Inversion. 

When  the  radius  of  the  crack  is  small  compared  to  the  wavelength  of  the  incident  TM  wave 

we  have,  for  small  £  =  . 

2 


/  _  n.  ,  \ 

B  e  /n! 

n 

[3m+  nj 

im2-  n2J 

l  2  J 

r~ 

ml 


(nc 


(49) 


and 


Ai  4 11  tl 


nil  —  sin 
(! 


<nc  *  (2iB_enm\\ 

>  -  I 


l2-n2 


(50) 


We  define  the  following 


m 

f  *  — —  sin  rtvp  nC.  (51) 

m  ml 


T 

mn 


3.m.t  n  , 

2  (m2-  n2) 


(52) 


(53) 


We  then  have 


* 

1 


T  b  -f  e 

mn  rt  m 


(54) 


(55) 
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Hence,  we  write 


A.  =  4  tl-— 
1  I! 


•  00  T  ^  f  E m 

I  T-  - 


n  m=l 


l2-  n2 


(56) 


where 


"  l2-  r>2 


For  small  f  (that  Is,  cracks  much  narrower  than  a  wavelength)  this  shows  that  the  only 
significant  At  is  Ay  ,  given  approximately  by 

2  ( nc 

»  -4  ne  9l  j  sin  <p  . 


We  see  from  the  above  that  the  diffracted  field  behaves  as  (ka)2  for  small  ka. 
* -Qn«(ka)2 H®\kr)  sin  ^sin  0<nc, 


(57) 


(58) 


and 


n  (ka)2H®  (kr)  sin  0sin  0<nc. 


(59) 


If  we  use  the  £ar-fleld  expressions  for  the  Hankel  functions,  we  write  for  plane  TM  wave 
scattering  from  a  small  cylindrical  crack  (  r  <*>). 


E<*I 


-9UV 2n  (ka)  e 


2  (3n/4  e 


-ikr 


Vkr 


sin  <p  sin  0 


Inc 


(60) 
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and 


(61) 


The  problem  of  plane  wave  scattering  from  a  slit  in  an  Infinitely  thin  conducting  plane 
has  been  solved,  and  Its  solution  Is  well  known.  It  has  been  discussed  in  a  recent  paper  14  and 
for  this  polarization  and  the  special  case  of  a  narrow  slit,  we  have 


J3n/4  e 


-Ikr 


ikr 


sin  <p  sin  <p 


Inc 


(62) 


which  Is  the  same  form  as  for  the  crack,  but  differing  by  a  constant  factor  — — — -  . 

4(^n) 

Performing  the  necessaiy  matrix  Inversion  we  find  that  =-0.19  which  gives  for  TM 
scattering  from  a  small  crack 


E^  =  0.19x(ka)2H|2)(ka)  sin  #sin  ^<nc. 


(63) 


To  check  the  consistency  of  the  small-crack  approximation,  we  have  written  a  computer 
program  to  solve  the  general-ka  Eqs.,  (45)  and  (46),  and  then  compared  the  general  results  to 
the  small-crack  results  as  ka  approaches  zero.  For  small  ka  the  first  mode  dominates,  and  a 
comparison  for  the  two  cases  is  shown  In  Figure  2.  It  can  be  seen  that  as  the  crack  becomes 
small,  the  closed-form  small-crack  results  coalesce  with  the  general  results. 


14  Shore,  R.A.,  and  Yaghjlan,  A.D.  (1988)  Incremental  diffraction  coefficients  for  planar 
surfaces,  IEEE  T  AP,  36(No.  1):55. 
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Figure  2.  General  and  Small-Crack  Q™  vs.  ka 


3.  TRANSVERSE  ELECTRIC  POLARIZATION 

3.1  Problem  Geometry 

We  consider  a  TE  plane  wave  that  is  polarized  in  the  z-direction  and  makes  an  angle  0 ,nc 
with  the  positive  x-semiaxis.  It  is  incident  upon  and  scatters  from  a  half-cylindrical 
indentation  in  the  ground  plane  (see  Figure  3).  The  crack  is  described  by  r  =  a  and  n  <  <f>  <  2n 
It  is  of  infinite  extent  in  the  z-direction  and  the  problem  is  two  dimensional.  All  angles  are 
measured  positive  In  the  counterclockwise  direction  starting  with  0=0  along  the  positive  x- 
semiaxis.  Angles  are  in  the  range  0  <  0  <2  k  . 
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Figure  3.  Problem  Geometry  for  Incident  TE  Waves. 

3.2  Exterior  Field 

The  incident  magnetic  field  is  a  plane  wave  written  as 


,lnc  (Icrcos  ($  (nc) 

*  ’  (64) 


where  unit  amplitude  has  been  assumed  and  e i0)t  time  dependence  suppressed.  Written  in 
cylindrical  Bessel  functions. 13  we  have 


£  <njn(ko 


(65) 


In  the  exterior  region,  the  scattered  field  will  be  considered  to  be  made  up  of  two  parts.  The 
first  is  the  reflected  wave  that  would  be  present  if  there  were  no  crack,  and  the  second  is  the 
deviation  from  this  caused  by  diffraction  at  the  crack.  Thus,  in  the  exterior  region  we  refer  to 
the  incident,  reflected,  and  diffracted  waves.  We  are  most  interested  in  the  diffracted  field. 

The  reflected  wave  is  a  plane  wave  and  is  well  known.  We  have  <f> re/  =  2n  -  <p inc  so  that  the 
reflected  plane  wave  can  be  written  as 


wr«f=  e(kr  cos  U+*<nC) 


(66) 
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In  terms  of  cylindrical  Bessel  functions  we  have 


i 


Cjjkr)  *  ** 


inc) 


(67) 


A  bit  of  algebra  and  a  few  trigonometric  identities  allow  us  to  show  the  sum  of  the  incident 
and  reflected  Hz- waves  satisfy  the  boundary  condition  of  zero  tangential  electric  field  on  the 
perfectly  conducting  ground  plane. 

For  the  diffracted  Hz- wave  we  write 


H^=  £  AnH^(kr)  cos  nip . 

rt*0 


(68) 


where  An  is  an  unknown  modal  coefficient,  and  we  have  expanded  this  outgoing  wave  in  terms 
of  cylindrical  Hankel  functions  (of  the  second  kind)  as  well  as  cos  n^.  Noting  that  because 
d#  cos  nip  =  -n  sin  nip  and  sin  0  =  sin  rur  =  0,  the  diffracted  Hz- wave  satisfies  the  boundary 
condition  on  the  ground  plane.  If  we  write  the  total  He-fleld  in  the  exterior  region  as 

H2=H<nc  +  Hre/+  Hd/,  (69) 

z  z  z  z 


it  is  clear  that  =  0  for  <p  =  0,  rr,  the  total  exterior  Hr-field  satisfies  the  boundary 

condition  on  the  perfectly  reflecting  ground  plane. 

We  now  proceed  to  write  the  incident  and  reflected  H*- fields  as  series  over  n  =  0  to  °®  rather 
than  over  all  n.  We  write  the  incident  field  as 


hT  =  ~J o(kr)+  i 

n=0 


ru»0 


(70) 


or 


H^c  =  ~JQ(kr)  +  2  X  <"Jn(kr)  cos  n(ip  -  <p <nc). 

nmO 


(71) 
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Similarly  for  the  reflected  field 


Hr*J  =  -JQ(kr)  +  2^  lnJn(kr)  cos  n  {<p  +  <p  <ne). 

n»0 


(72) 


Hence  we  can  write,  using  the  ldentiy  cos  n 


(# 


(nc\  (nc  (nc 

ip  )=  cos  nip  cos  nip  t  sin  n#  sin  nip 


+  H^=  -2JQ(kr)  +  4  £  inJn(kr)  cos  nip  cos  nip  <n°. 

rt=0 


(73) 


Recalling  our  expression  for  the  diffracted  field,  we  write  the  toted  exterior  Hz-field  as 


H2  =  -2JQ(kr)  +  X  [4  <n^n(^r)  cos  >nC  +  An/f|2)(fcr)|  cos  nip . 

nsO 


(74) 


3.3  Interior  Field 


In  the  Interior  region  the  total  Hz- field  can  be  written  as 


(75) 


where  Dn  is  an  unknown  modal  coefficient.  However,  for  this  calculation.  It  will  be  more 
convenient  if  we  write  this  equivalently  as 


hJ  =  X  •dn(>cr)  (Bn  cos  nip  +  Cn  sin  nip), 

rtm  0 


(76) 


where  B„  ,  Cn  are  the  modal  coefficients  that  will  be  determined  from  the  boundary  conditions. 
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3.4  Electric  Field  Component 

The  ^-component  of  the  electric  field  Is  also  of  Interest  In  this  problem.  From  Maxwell's 
equations  we  have 


=  -i-  d  H 

cue  r  *  (77) 


so  that 


E2  =  Jk- 


coe 


{- 


2  J'0(kr)  + 


nmO 


4  fn  J'(kr)  cos  n<t><nC  +  A  H®  (kr) 

rv  '  n  n  v 


cos  n0 


} 


and,  for  clarity 


Y  A  H{2V(kr)  cos  n«J. 
0  .3,  rl  n 


3.5  Boundary  Conditions 

Boundary  conditions  required  for  solving  this  TE  problem  are 


E^1  =  0  for  r  =  a  and  n  <  <p  <  2n , 
1  2 

Eq  =  Ep  for  r  =  a  and  0  <  <p  <  n , 
1  2 

H  -H  for  r  =  a  and  0  <  <p  <  n . 

Z  Z 


(78) 


(79) 


(80) 


(81) 

(82) 

(83) 
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The  first  condition  enforces  zero  tangential  electric  field  on  the  surface  of  the  crack  and  the 
last  two  ensure  continuous  tangential  electric  and  magnetic  fields  across  the  aperture.  Recall 
that  the  total  exterior  field,  as  written,  satisfies  the  boundary  conditions  along  the  ground 
plane. 

Applying  the  boundary  conditions,  we  get  the  following  three  equations 


X  BnJ  C0S  n^l  +  X  Cn  J'n(ka)  sln  n(t>i  =  °. 

n=0  rtxO 


(84) 


X  BnJn(ka)  COS  n<t>2+  X  C,  Sin  TVp2 

rb=0  ns=0 


=  -2  J'0(ka)  +  X  4  cos  n<j> 


n=0 


inc  +  AnH®]  (ka)j cos  n<p2. 


(85) 


X  B^J^ka)  cos  n<?2  +  X  CnJn(ka)  sin  n02  (86) 

rt*0  n=0 


=  -2  JQ(ka)  +  X  I4  cos  ^  +  AnHj21(ka)]  cos  n<p2 

n=0 


where  /r<  <  2/r  and  0  <  4>2  <  n  .  We  now  define  the  following  for  convenience 


Fn  =  {(4  <"  cos  ntf '"c -  Bj  Jn(ka)  +  ^H^ka) }, 


(87) 


G  =-B  J  (ka), 
n  n  n 


(88) 


R  =  C  J  (ka) 

n  n  n 


(89) 


as  well  as 


dFn 

d(ka) 


=  F'n  and  so  on.  The  three  Boundary  Condition  equations  thus  become 


X  (-  l)n  G'n  cos  =  X  (“  l)nB  'n  sin  n<p  +  BQ  J'0(ka),  (90) 

rm  1  na  1 


21 


where  we  have  made  the  change  of  variables  4>x  ->  0  +  n  and  02  -»  0 . 


X  F'ncosn<p  =  '£  R'nsinn0+(2-Bo)J'o(ka)-AoH^)(ka).  (91) 

n*l  r>=l 


X  F  cos  n0  =  X  R  sin  n0  +  (2  -  BQ)  JQ(ka)  -  A()f^t2)(ka) .  (92) 

rvx  1  n=  1 


All  three  equations  are  valid  for  0  <  0  <  n  and  all  are  of  the  same  form  with  respect  to  0 . 


3.6  Orthogonality 

We  consider  the  equation 


£  Fn  cos  nip  =  £  sin  nip  +  K(ka) .  (93) 

rt=  1  n*l 


where  Fn  .  Rn  ,  K(ka)  are  not  functions  of  0  and  0  <  <p  <  rr  .  We  have  the  following  orthogonality 
relations  among  the  sinusoids 


cos  nip  cos  m0  dtp  =  8 

2  mn 


(m,n>  1), 


(94) 


and 


sin  n0  cos  mtp  dip  = 


3  B  -  if  n  -  m  is  odd 
_2  2 
n  -  m 

0  if  n  -  m  is  even 


(95) 


We  also  note  that 
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cos  m<p  d<p  =  0 


m  =  1,  2,  3  ... 


(96) 


Now,  we  multiply  Eq.  (93)  by  cos  m^and  Integrate  from  0  t on. 


r 


rtm  l 


Fn  cos  rup  cos  m<p  d<p 


mK  OO 

IK 

Jo  n»l 


sin  nip  cos  m<p  dip  + 


f 


cos  mipdip. 


(97) 


The  above  relations  allow  us  to  write 


F 

m 


I 


(  ) 


(m£  1). 


(98) 


as  a  solution  to  Eq.  (93).  Here  X  n  Indicates  sum  over  n  from  1  to  ~  for  those  values  of  n 
where  (n  -  m)  is  odd. 

The  case  m  =  0  is  considered  separately.  In  this  case  we  effectively  multiply  Eq.  (93)  by 
unity  and  then  Integrate  from  0  to  n.  We  note 


sin  nip  dip  = 


Z 

n 

0 


n  odd 
n  even 


(99) 


We  then  have.  In  addition  to  Eq.  (98)  the  relation 


M 

X 

rw  1.3.5 


I  r  =-*K  (m  =  0). 
n  n  2 


(100) 


The  three  Boundary  Condition  equations  thus  have  the  solutions,  for  mS  1, 
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G' 

m 


=  4X 


L(  nR' 


2  2  . 
n  -  m  J 


(101) 


F' 

m 


I 


(  "*«  1 


(102) 


f  n« 


\ 

> 


(103) 


or,  recalling  our  previous  definitions 


B  J'  (ka) 

mm7 


(nC  J'(ka)) 

n  rr  ' 

2  2 
V  n  -  m 


( 104 ) 


(4  r  cos  rrup  flj  J'm(ka)  +  AmH®  (fca )  -J2 


nC  J'(ka) 


2  2 

n  -  m 


(105) 


(4  tm  cos  m<f>  B J  Jjka)  +  AJ^(ka)  =  ££ 


nC  J  (ka) 


2  2 

n  -  m 


(106) 


For  m  =  0  we  have 

a* 

y  i-C  J'(ka)  =  S8  J'(ka), 

rar3.5 n  n  n  2  0  0 
f  HCnJ'n<fca)=f  {(2"B0)  ^fca)+A0H02,  (fca)}. 
£  i  c„  •'<'“>>  *  ff  { (2-Bo)  ♦  AflSV0)}  ■ 

n*  1,3.5  z  L  J 


(107) 


(108) 


(109) 


These  two  sets  of  equations  will  be  used  to  solve  for  the  modal  coefficients  An  ,  Bn  ,  and  C„  . 
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3.7  Solution  for  Modal  Coefficients 


We  first  eliminate  Am  and  A0  respectively  from  the  two  sets  of  equations  above,  and  then 
we  use  the  Wronskian  relationship,  W|ym(ka)H^(ka)]  = — ,  to  write 


, 


J'„(fca)Hm!(ka)  ~  JrSka)H™'(ka) 


+  -2J_ B  -8-(-L\ 

n  ka  m  n  \ka / 


m  Inc 

I  cos  m<t>  , 


(110) 


where  we  have  now  combined  the  two  sets  of  equations.  Substituting  for  Bm  ,  and  rearranging, 
we  write,  for  m  =  0,  1.2,  3... 


n  im  J'n  (ka)  cos  mtp >nc. 


which  is  an  expression  that  we  can  solve  for  Cn  .  Then,  when  we  have  Cn  we  can  write 


(112) 


and,  for  l  =  1,  2,  3  ... 


A 


l 


8  y 

( nC  J'  (ka) 

n  rr  ' 

[*i 

2  .2 

n  - 1  / 

-4  ll  J’lka )  cos  10^ 


ffy(ka) 


(113) 


3.8  Small  Crack  Limit 

When  ka-*  0  the  cylinder  functions  have  the  following  simple  forms 
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J((fca)  -> 


(ka/  2)1 
(! 


(114) 


and 


ff\ka) 


J _ IL_ 

(ka/2)1 


1*0;  H*2)(ka)  ->  -2-  ln(ka). 

°  iff 


(115) 


The  equations  that  we  must  solve  also  become  much  simpler,  and  we  can  separate  the  ka- 
dependence  from  the  necessary  numerical  matrix  Inversion. 

When  the  radius  of  the  crack  Is  small  compared  to  the  wavelength  of  the  Incident  TE  wave 
we  have,  setting  m  =  0  In  Eq.  (Ill)  and  then  substituting  Into  Eq.  (1 12), 


A  =  iS-(ka)2 
o  2 


1  +  i  (ka)2  ln(ka) 
2 


(116) 


and,  with  e  =  ka/2  «  1  and  c  =C  J'(ka), 

n  n  n 


auii 


*  (  2nc  \ 
1 


vn2-  l2 


-  £  il  -L  —  cos 
2  He 


(nc 


0*1). 


(117) 


where  the  cn  are  found  from 


y  ! 3.11 1  in']  c  *  ff  cos  m<p,nc  (m,  rt  2  1) . 

nin2-m2J  n  m!  £ 


If  we  define 


=  ff  ( 


m 


JB-  cos  m<p <nc, 
ml 


(119) 


T 

mn 


3ii±.fn 
2  2 
n  -  m 


(120) 
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we  have 


V  T  c  =/  l  (m,n>l), 

mn  n  ~'m  v  7 


Hence,  we  write 


l  *  •  a  —,rr~  i  f  i  • 

II  f-f  'If  =4<7Fl9im-t 

*•  n  m  ri  —  l  1  m 


where 


_  .  4  nT  „ 

Q  =  V  _ QSl-  s  . 

9lm  **  n  2  ,2  lm 

n  -  I 


For  small  e  (that  is.  cracks  much  narrower  than  a  wavelength)  this  shows  that  the  only 
significant  A,  are  A0  ,  A;  .  with  Ax  given  approximately  by 


2  _  Inc  „  . ,  .2  ,  fnc 

Aj  *-4j£  cos  <t>  =-9ux(ka)  cos  <p 


We  see  from  the  above  that  the  diffracted  field  behaves  as  (ka)2  for  small  ka. 


Hdif~-£(ka)2  <f-l+ - 2- -  H®(kr)  +  2Q  H®\kr)  cos  ^cos  .  (124) 

*  2  1  +l(ka)2  In  (ka) 

L  V  2  J 


If  we  use  the  far-field  expressions  for  the  Hankel  functions,  we  write  for  plane  TE  wave 
scattering  from  a  small  cylindrical  crack  (r  -»  «>) 


—  (icr  p 

-  -  -i/ff"  ( ka)2  ei3n/*  e-—r  -1  + - - +  20  cos  0  cos  <p 

‘  '  2  '  Vk?  1  +  i-  (ka)2  In  (ka)  M 
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where 


(126) 


Performing  the  necessary  matrix  Inversion  we  find  that  Q\  i  =  0.18  which  gives  for  TE  far-fleld 
scattering  from  a  small  crack 


H^f 


{2rt  (ka)2 


13W 


fkr 


1  +  i  (ka)2  In  (ka) 


(127) 


A  computer  program  has  been  written  to  check  the  consistency  of  the  small  crack 
approximation.  The  general  expressions.  (Eqs.  (Ill)  -  (113)),  are  solved  numerically  and  then 
compared  with  the  small-crack  results.  Figure  4  shows  the  negative  of  the  bracketed  quantities 
from  Eq.  (1 16)  compared  to  the  general-ka  expressions  as  ka  becomes  small.  Note  that  for  ka 
very  small,  the  bracketed  terms  in  Eq.  (116)  tend  to  negative  unity.  Figure  5  shows  the 
coefficient  Qi  \  compared  to  the  lowest  mode  coefficient  from  the  general  program  as  ka 
becomes  small.  For  small  values  of  ka,  (where  (ka)2  can  be  neglected  compared  to  unityl  the 
lowest  mode  dominates,  and  the  general  and  small-crack  coefficients  of  this  mode  agree. 

The  problem  of  plane  wave  scattering  from  a  slit  In  an  infinitely  thin  conducting  plane 
has  been  solved  and  is  well  known.  For  the  TM  polarization,  we  found  the  results  for 
scattering  from  the  small  semicircular  crack  to  be  in  agreement  with  those  for  the  silt 
(differing  by  a  constant  factor.)  For  TE  scattering  from  a  narrow  slit.  Reference  [14]  gives 


H*'-  V  27r  e <3n/4 


e~ikr  1  +  .25(ka)2(sin2  finc  -  cos2  pi) 
Vkr  2  In  (ka/4)  +  2y  +  tn 


+  0.25  (ka)2 


(nc 

COS  p  COS  p 


(128) 


where  yis  Euler’s  constant.  We  note  that  for  the  TE  polarization,  the  leading  term  for  the 
small  crack  is  (ka)2  while  the  leading  term  for  the  slit  is  jn  we  conclude  that  for  the 

TE  polarization  the  narrow  slit  scatters  more  than  the  small  crack.  The  physical  reasoning  is 
as  follows.  For  TM  polarization,  the  currents  excited  by  the  Incident  electric  field  are  normal 
to  the  crack/sllt.  In  the  former  case  the  small  crack  and  the  small  slit  appear  to  the  excited 
currents  to  be  similar,  while  in  the  latter  case  the  currents  are  stopped  completely  by  the  open 
circuit  of  the  slit,  but  not  by  the  perfectly  conducting  crack. 
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KA 

Figure  5.  General  and  Small-Crack  Q vs.  ka. 

4.  SUMMARY 


In  this  report  we  have  developed  expressions  for  the  fields  scattered  from  a  semicylindrical 
crack  in  a  perfectly  conducting  ground  plane  for  both  transverse  magnetic  and  transverse 
electric  polarizations.  Although  a  numerical  matrix  inversion  is  required  in  the  solution, 
when  the  crack  is  much  narrower  than  the  wavelength  of  the  incident  radiation,  the  ka- 
dependence  is  removed  from  the  numerical  matrix  inversion.  The  closed-form  results  that  are 
found  for  small  cracks  are  thus  functions  of  ka  and  are  valid  for  any  value  of  ka  such  that  ka 2 
can  be  neglected  compared  to  unity. 

Since  the  TM  and  TE  results  for  small  cracks  are  functionally  dependent  upon  ka.  they  can 
be  combined  vectorially  to  consider  off-normal  scattering  from  small  cracks.  Because  of  this, 
the  results  presented  here  for  scattering  from  a  small  straight  crack  can  be  extended  by  such 
methods  as  Incremental  Diffraction  Coefficients  (IDC)  14  to  calculate  the  scattering  from  a 
small  curved  crack.  Large  reflector  antennas  are  typically  assembled  from  several  sections, 
leaving  gaps  between  them.  Aircraft  have  gaps  around  doors  and  access  panels.  The  former 
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alfect  the  radiation  pattern  of  the  antenna  and  the  latter  contribute  to  the  aircraft's  radar 
cross  section.  Using  IDC,  Shore  and  Yaghjian  14  have  modeled  these  as  slits  in  thin 
conducting  screens,  but  the  geometry  of  a  cylindrical  crack  seems  a  more  physically  realistic 
model.  Since  the  crack  and  slit  only  agree  well  for  the  (small  ka)  TM  polarization,  real  cracks 
may  be  better  approximated  using  the  results  presented  here. 
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MISSION 

°f 

Rome  Air  Development  Center 


RADC  plans  and  executes  research,  development,  test  and 
selected  acquisition  programs  in  support  of  Command,  Control, 
Communications  and  Intelligence  (CSI)  activities.  Technical  and 
engineering  support  within  areas  of  competence  is  provided  to 
ESD  Program  Offices  (POs)  and  other  ESD  elements  to 
perform  effective  acquisition  of  C3I  systems.  The  areas  of 
technical  competence  include  communications,  command  and 
control,  battle  management  information  processing,  surveillance 
sensors,  intelligence  data  collection  and  handling,  solid  state 
sciences,  electromagnetics,  and  propagation,  and  electronic 
reliability /maintainability  and  compatibility. 


